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Abstract 



Context. In recent years new observations of pre-main sequence stars (pre-MS) with Z < Z© have been made available. To take full 
advantage of the continuously growing amount of data of pre-MS stars in different environments, we need to develop updated pre-MS 
models for a wide range of metallicity to assign reliable ages and masses to the observed stars. 

Aims. We present updated evolutionary pre-MS models and isochrones for a fine grid of mass, age, metallicity, and helium values. 
Methods. We use a standard and well-tested stellar evolutionary code (i.e. FRANEC), that adopts outer boundary conditions from 
detailed and realistic atmosphere models. In this code, we incorporate additional improvements to the physical inputs related to the 
equation of state and the low temperature radiative opacities essential to computing low-mass stellar models. 

Results. We make available via internet a large database of pre-MS tracks and isochrones for a wide range of chemical compositions 
(Z = 0.0002 - 0.03), masses (M = 0.2 - 7.0 Mq), and ages (1 - 100 Myr) for a solar-calibrated mixing length parameter a (i.e. 1.68). 
For each chemical composition, additional models were computed with two different mixing length values, namely a = 1.2 and 1.9. 
Moreover, for Z > 0.008, we also provided models with two different initial deuterium abundances. The characteristics of the models 
have been discussed in detail and compared with other work in the literature. The main uncertainties affecting theoretical predictions 
have been critically discussed. Comparisons with selected data indicate that there is close agreement between theory and observation. 

Key words. Stars: pre-main sequence - Stars: evolution - Stars: formation - Stars: interiors - Stars: low-mass ~ Stars: HR diagrams 



1. Introduction 

We present a new set of pre-main sequence (pre-MS) stellar 
tracks and isochrones based on both the state-of-the-art input 
physics and outer boundary conditions. These models cover a 
large range of metallicities, helium abundances, masses, and 
ages, providing useful tools for interpreting observational data. 
Pre-MS tracks and isochrones collectively represent the theoret- 
ical tool needed to infer the star formation history and the ini- 
tial mass function of young stellar systems. Thus, the availab- 
ility of a large database with a fine grid of chemical compos- 
itions, masses, and ages will allow us to improve our under- 
standing of the star formation process taking full advantage of 
the growing amount of data o f very young clusters and associ- 
ations in the Milky Way (e .g. IStolte et al.ll2005t iDel^ado et aP 



2007"; 'Brandner et al."2008'), the Small Magellanic Cloud (e 
Gouliermis et al. 2006b; Nota et al. 2006; C arlson et al. 2007|; 
Gouliermis et al. 2007b; Sabbi et al.l 120071; iGoulierm is et al.l 
2008.; iCignoni et al. 2 009.. ,20 id. and the Large Magellanic 
Cloud (e.g.lR omanie Uo et al.1 120041: iGoiil iermis et al. 2006i 
Romaniello et all l2006nGouliermis et all l2007at IDa Rio et'all 
20091) ! 



Send offprint requests to: P. G. Prada Moroni, prada@df.unipi.it 

* Tracks and isochrones are available on the web at the URL: 
|http ://as tro.df.unipi.it/stellar-models/ 

** Tracks and isochrones are also available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+A/= 



As a n example of a n intriguing application of pre-MS iso- 
chrones, ICignoni et al.l (2009) inferred the star formation his- 
tory of NGC602, a very young stellar system in the Small 



Magellanic Cloud (see also'Nota et al.'2006l: ICarlson et al.l2"007t 
Gouliermis et al. 2007b; Sabbi et al. 2007). 

From the th eor etical point of view, as ear ly understood by 
Hayashi" ('1961') and'Havashi & Nakanol (ll963h . the evolution of 
hydrostatic pre-MS stars is in principle less complex than more 
advanced evolutionary phases from the numerical point of view, 
since no special algorithms are required to follow what it is es- 
sentially a quasi-static contraction. However, the pre-MS track 
location in the HR diagram strongly depends on the physical 
ingredients used in the evolutionary codes, such as the equa- 
tion of state (EOS), the low-temperature radiative opacity, the 
outer boundary conditions , and the adopted convection treat- 
ment (see e.g.. | D'Antonaj[r 993; D'Antona & Mazzitelli 199|, 
'1997; Baraff^e et al. 1998; Siess et al. 2000; Baraff^e et al. 20ol 
iMontalban et al..2004.) . The uncertainties due to these quantities, 
progressively increase as the stellar mass decreases, especially 
for very low-mass stars (i.e. M $ 0.5 Mq). Therefore, to com- 
pute models as accurately as possible, it is mandatory to include 
the most recent updates of the above quoted physical inputs. 
In the past two decades, several generations of theoretical pre- 
MS models have been developed following improvements to the 
physical inputs to provide progressively more accurate theoret- 
ical pr edictions (e.g. Baraffe et al. 1998; D' Antona & Mazzitelli 
1997] 
120091) 



ISiess et al.1 [20001: iDotter et al.1 120071: fPi Criscienzo et al 
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The theoretical tracks and isochrones for pre-MS stars de- 
scribed in the present paper have been computed with an updated 
version of the FRANEC e volutionary code (Degl'Innocentiet al] 
l2008t IValle et aH l2009h relying on the most recent phys- 
ical inputs. The models cover a wide range of metallicity 
(Z = 0.0002 - 0.03), masses (M = 0.2 - 7.0Mo) and ages 
(t = 1 - 100 Myr). The corresponding database is available on 
the wetQ. 

The input physics adopted in the calculations are described 
in Sect. 2, followed by a discussion of the main sources of uncer- 
tainty in the evolution of young stars. In Sects. 3 and 4, our tracks 
and isochrones are compared with recent evolutionary models 
and data. Finally, in our last section the database is described. 



2. The models 

We describe the adopted physical inputs focusing only on the 
ingredients that affect mainly the pre-MS evolution and the loc- 
ation of the track in the HR diagram. 

Before starting this description, we recall that our mod- 
els are standard in the sense that they do not include rota- 
tion, magnetic fields, and accretion. This is an important point 
to emphasize because the inclusion of these processes affects 
the evolution of pre-MS models introducing additional sources 
of uncertainty (s ee e.g., |Pinsonneault et al. 1990; Siess et al., 



^1997; D'Antona erai]|200d; Siess 2001; Stabler & Pallal 120041: 
Baraffe et al...2009HHosokawa et al...201 1.) . 

Although inconsequential to pre-MS evolution, we also note 
that our current evolutionary tracks take into account micro- 
scopic diffusion. In particular, the diffusion velocities of He, Li, 
Be, B, C, N, O, and Fe are directly co mputed following the 
method described in iThoul et al] d 19941) . while the other ele- 
ments are assumed to diffuse as Fe but with a velocity scaled 
to their abundances. We do not consider radiative levitation. 



2.1. Equation oi State 

To compute stellar models, it is necessary to use an accurate 
EOS covering the very wide range of temperatures and dens- 
ities spanned during the entire evolution. We adopted the most 
updat ed version of the OPAL EOSg in the versi on released in 
2006 ([Rogers et al.|[T996l;lRogers & Navfonovl2002i) . 

Figure [T] shows the validity domain in the temperature- 
density plane of the OPAL EOS. Thick solid lines represent the 
time evolution during the pre-MS phase of the stellar centre (top- 
right side) and the atmosphere base at an optical deptlQ rph=10 
(bottom-left side) of the 0.2 Mq and 7 Mq models with Z - 0.03, 
Y = 0.290, and a - 1.2, which is our coldest set of tracks. To 
clearly show that the OPAL EOS is fully capable of covering 
the space of parameters required to compute all the models in- 
cluded in the database, we added to Fig.[T]the time evolution of 
the central and outer regions of our hottest models, that is, those 
with Z = 0.0002, Y = 0.250, and a = 1.9 (thick dashed Unes); 
in this case, the entire structure profiles for the starting model at 
the beginning of the pre-MS evolution (dotted lines) and for the 
final ZAMS model (dot-dashed lines) are also shown. 



' Models and isochrones are available in the theoretical plane and 
will soon be available in different photometric systems. 

2 The OPAL EOS tables are available at the URL: 
[FEt p : //opalopacity . llnl . gov/EOS_2885/ 

^ As we will describe in Sect. l2.2l we choose Tpi, = 10 for the match- 
ing point between the atmosphere model and the interior. 
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Figure 1. Validity domain of the OPAL EOS (shaded area) in 
the plane (log T, log p); we superimposed the time evolution of 
the bottom of the atmosphere (xph^lO) and the stellar centre for 
Z = 0.03 (thick solid line) and Z — 0.0002 (thick dashed line), 
and the whole structure at the beginning of the contraction along 
the Hayashi track (dotted line) and in the ZAMS (dot-dashed 
line) forZ^ 0.0002. 
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Figure 2. HR diagram for evolutionary tracks in the mass range 
0.2 - 2.0 Mo computed adopting the OPAL EOS 2006 (solid 
line) and the FreeEOS in the PTEH configuration (dashed line) 
for the labelled chemical composition and a. 



Figure |2] shows the effect of two different EOS widely used 
in the literature, nar nely the OPA L EOS 2006 and the FreeEOS 
(for more details see llrwinl2004l) in the confi guration that should 
reproduce the PTEH EOS (i Pols et al.|[T995i) . The effect of the 
EOS has been widely disc ussed in the literature in the case 
of pre-MS objects (see e.g, iMazzitelhl [19891 iD'Antonal [19931 
and refer ences therein) and for low-mass main sequence stars 
(see e .g., lDormanetal.lll989t lNeecelll984l: IChabrier & Baraff^d 
1 1997'); here we limit our discussion to presenting the results of 
our comparisons. For stars more massive than 1 M©, the effect is 
essentially limited to the Hayashi track; the difference in effect- 
ive temperature is about 50 K for 1 M© model and progressively 
decreases as the mass increases becoming about 30 K for 2 Mq 
and almost negligible for more massive stars. For models less 
massive than 1 Mq, the EOS also affected the ZAMS because 
of thick convective envelopes present in the hydrogen burning 
phase, whose structure depends on the adiabatic gradient. For 
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Table 1. Differences between models computed adopting the 
PTEH and the OPAL EOS. The second column shows for each 
value of the mass the maximum difference at a fixed logL/Lg 
along the Hayashi track. The third and forth columns show the 
differences in Tgff and log L/Lq, respectively, for ZAMS models. 



M/Mo 


AT 


eff 


AlogL/Lo 




Hayashi 


ZAMS 


ZAMS 


0.2 


-30 K 


-30 K 


-0.02 


0.4 


-35 K 


-48 K 


-0.04 


0.6 


-40 K 


-60 K 


-0.03 


0.8 


-46 K 


-35 K 


-0.01 


1.0 


-48 K 


-15 K 


- 


1.5 


-34 K 


- 


- 


2.0 


-30 K 


- 


- 



this mass range, the maximum differences between the tracks 
with OPAL and the PTEH EOS at fixed luminosity, shown in 
Table [U are of about 30 - 40 K along the Hayashi track and 
30 - 60 K in ZAMS depending on the mass. We note that for 
low-mass stars the ZAMS luminosity also depends on the adop- 
ted EOS. 



2.2. Boundary conditions 

To solve the differential equations of the stellar interior structure, 
suitable outer boundary conditions are required. The usual ap- 
proach followed in standard evolutionary codes consists of tak- 
ing as outer boundary conditions the values of the main phys- 
ical quantities at the base of the atmosphere. These quantities 
are provided by a direct integration of the equations describ- 
ing a mono-dimensional atmosphere both in hydrostatic equi- 
librium and for the diffusive approximation of radiative trans- 
port, in addition to a grey T(t) relationship between the tem- 
perature and the optical depth r. The T(t) relationships most 
commonly chosen are t hose of the Eddington approximation and 
iKrishna Swamvl (1 19661 hereafter KS66). 

_However, as discussed in several papers (see e.g., Auman 
1969; Dormanetal. 1989; Saumonet al. 1994; B araffe et al 
1995; Allard et al. 1997; Ch abrier & B araffe 1997; Baraffe et all 
1998t i200Z) . this technique is too crude for low-mass stars and 



in general for cold atmospheres. A much more sophisticated ap- 
proach consists of adopting as outer boundary conditions for 
the integration of the stellar interior equations the main physical 
quantities at a given optical depth provided by detailed, non-grey 
atmospheric models that solve the full radiative transport equa- 
tion. 

We adopted the atmospheric models by Brott & Hauschildj 

2005 h. herea fter BH05, c omputed with the PHOENIX code 

Hauschildt & Baronl|1999h . that are available for the parameter 

ranges 3000 K < Teff < 10 000 K, -0.5 < logg(cms-2) < h-5.5, 

and -4.0 < [M/H1 < -hO.5, adoptin g aatm - 2.0 and the solar 

mixture of iGrevesse & Noelsl (Il993h . 

In the range 1 0000 K < Teif < 50000 K, 
4-0.0 < logg(cms-2) < H-5.0, and -2.5 < [M/H] < H-0.5, 
where the models of BH05 are unavailable, we used the 
atmospheric structures of iCastelli & Kurucz (2003), hereafter 
CK03. These models are computed for the solar mixture of 
IGrevesse & Sauva]|(ll998h . As in the previous case, the mixing- 
length scheme is followed to describe the convection, but with 
a lower value of the mixing-length parameter aatm = 1.25. The 
effect of the adoption of two different o-atm values should be 
small compared to the other uncertainty sources, since for ef- 
fective temperatures above 10000 K, which roughly correspond 



to ZAMS stars of M <; 2.0 M©, the convective envelope, when 
present, is very thin. 

From these two sets of atmospheric models, we built tables 
of boundary conditions for a fixed value of Tph, defined as the 
optical depth at which the interior and atmosphere matches. In 
more detail, these tables contain the temperature T(Tph, Tetf, 
log g, Z) and the pressure P(Tph, Teff, log g, Z) at the optical depth 
Tph extracted from the atmospheric models for the available Tgff, 
log g, and Z. During the computation of the stellar structure and 
evolution, a spline interpolation of these tabulated values at the 
Teff, log g and log Z of the star is performed to obtain the corres- 
ponding T(rph) and P(Tph) required for the interior integration. 
For our reference set of models we chose the commonly sugges- 
ted value of Tph = 10 jMorel et al.lll994t) ; this value represents a 
good compromise between two opposite requirements, namely, 
a large Tph to guarantee the validity of the photon diffusion ap- 
proximation, and a low Tph to avoid large discrepancies intro- 
duced by the inconsistencies often present between atmospheric 
and interior models. We note that we performed the interpolation 
in the global metallicity Z without taking into account the differ- 
ent heavy element distributions in the interior and atmospheric 
regions. 

It is worth noting that the interpolation technique adopted 
to obtain the boundary conditions must be chosen carefully, be- 
cause the available atmospheric models are not usually com- 
puted with this in mind and the resolution in Teff, logg, and 
[M/H] is often not very high. This might cause problems, mainly 
in the region of the HR diagram populated by red giant stars, in 
which the thickness of the convective outer layer changes ab- 
ruptly varying Teff and log g. As an example, we verified that if a 
spline interpolation is replaced with a linear one, large and non- 
physical waving appears in the HR diagram along the Hayashi 
phase. 

2.3. Opacity 

We adopted th e version of the OPA Ip radiative opacity 
tables (see e.g. Ilglesias & Rogerslll996h released in 2005 for 
log T(K) > 4.5; for this range of temperatures, the radiative 
opacity and the EOS are fully consistent, since they have been 
computed by the same group adopting the same physical pre- 
scriptions. For lower temperatures, characteristic of the outer 
stellar regions, we used the radiative opacity of Ferguso n et al] 
(l2005h (hereafter F05Jj, which include an accurate description 
of the molecular absorptions; indeed for T < 4500 K, hence for 
low-mass stars, molecules become the main source of opacity. 
We also mention, for the sake of complete ness, that we used the 
conductive opacity of Potekhin (1999) and Shternin & YakovlevI 
(12006) (hereafter PS Y06o however, for the present calculations 
thermal conduction is completely negligible, and we verified that 
even our lowest mass models are unaffected by increasing the 
conductive opacity of 5%. 

Our reference set of model s was computed adopt ing the solar 
heavy element distribution of lAsplund et al.l (l2005l) in both the 
low and high temperature regimes. 

Figure [3] shows the range of validity in the temperature- 
density plane of the F05, OPAL, and PSY06 opacity tables, we 



The OPAL radiative opacity can be found at the URL: 



http : // opalopacity . llnl . gov/new . html 

' The F05 low temperature radiative opacities are available at the 
URL: http : //webs . wichita . edu/physics/opacity 

* PSY06 conductive opacities are available at the URL: 
http : //w ww ■ iof f e . ru/astro/conduct/ 1 
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Figure 3. Validity domain of the radiative opacity tables (OPAL 
and F05, shaded area) and of the conductive ones {solid box). 
As in Fig. [T] we superimposed the time evolution of the bottom 
of the atmosphere (Tph=10) and the stellar centre for Z = 0.03 
(thick solid line) and Z = 0.0002 {thick dashed line), and the 
whole structure for the first model on the Hayashi track {dotted 
line) and for the ZAMS {dot-dashed line) for Z = 0.0002. 



note that the actual validity domain of the OPAL tables extends 
down to logT(K) = 3.75, partially overlapping with the F05 
tables, but we do not show it in figure, since we used the F05 opa- 
cities for log T(K) < 4.5. As in Fig.[Tl we plotted the time evolu- 
tion during the pre-MS phase of the stellar centre (top-right side) 
and the atmosphere base at an optical depth Tph=10 (bottom-left 
side) of the 0.2 Mg and 7 Mq models with Z = 0.03, Y = 0.290, 
and a - 1.2, i.e. our coldest set (thick solid line). The figure also 
shows the time evolution of the central and outer regions of our 
hottest models, that is, those with Z = 0.0002, Y - 0.250, and 
a - 1.9 (thick dashed lines); in this case, the entire structure pro- 
files for the starting model at the beginning of the pre-MS evol- 
ution (dotted lines) and for the final ZAMS model (dot-dashed 
lines) are also shown. 

We note that the current version of the radiative opacity 
tables do not extend sufficiently into the temperature-density 
plane to cover the entire pre-MS evolution of stars less massive 
than 0.6 - 0.5 M© (the exact value depending on the chemical 
composition), hence extrapolation to higher densities is required 
for both OPAL and F05 opacity tables. To our present know- 
ledge, there are no opacity tables available in the literature that 
cover the whole range of temperatures and densities spanned by 
low-mass pre-MS models. Moreover, the extrapolation of radi- 
ative opacities is a dangerous procedure, because opacity coef- 
ficients are very sensitive to temperature and density. Given this 
situation, it is of primary importance to evaluate the sensitivity 
of the computed evolution to the adopted extrapolation method. 

We tested three diff'erent methods: constant, linear for the last 
two points, and a linear fit to the last four points extrapolation. 
The reason for trying this latter method is to extend the mean 
slope of the opacity with respect to the density to reduce the 
contribution of a possible jump in the opacity near the border. 

We performed test of both non-grey and grey models, the 
latter ones ha ving been compu ted adopting the classical T(t) re- 
lationship of [Krishna Swamvl ([1966). We found that the extra- 
polation method does not aff'ect the non-grey models, because 
the extrapolated values of the radiative opacity are used only in 
a very thin shell, but not for computing the boundary conditions, 
which are provided by atmospheric models. In contrast, low- 



1 - J.5M 







Two-point linear extrapolation 
Four-point linear extrapolation 
Constant extrapolation 



Grey models, Y=0.278, Z=0.015, a=1.68 



3,8 



logT^ff(K) 



Figure 4. Comparison in the HR diagram of low-mass tracks 
computed adopting three different methods of opacity extrapol- 
ation (see text) in the regions not currently covered by opacity 
tables (low temperature-high density, see Fig.[3]|. 



mass grey models actually use the extrapolated radiative opa- 
city coefficients over the whole outermost region, atmosphere 
included, making the structures much more sensitive to the ad- 
opted extrapolation method. 

Figure |4] shows the HR diagram of grey pre-MS models 
computed with the three different extrapolation techniques. The 
impact of the extrapolation is completely negligible for stars 
more massive than 0.6 M©, because their structures are entirely 
covered by the OPAL+F05 tables, and becomes progressively 
larger and larger as the mass decreases, becoming quite import- 
ant for stars less massive than 0.4 Mq. The constant extrapola- 
tion, which is the most crude possibility, leads to ZAMS mod- 
els hotter than those computed with a linear fit extrapolation of 
about 170 K and 160 K for 0.4 Mq and 0.2 Mq, respectively. The 
difference between linear extrapolation with two or four points 
is lower (by about 45 - 60 K). 

2.4. Solar mixture of heavy elements 

A quite important, although often not explicitly specified, in- 
gredient of stellar modelling is the heavy element mixture, which 
is the number distribution of chemical elements heavier than he- 
lium. The adopted mixture affects both the radiative opacity of 
the stellar matter and the nuclear burning efficie ncy, through, for 

example, the C, N, O abundances (see e.g. Deg lTnnocenti et all 

II 1 1 ~ 

2006; Se stito et al. ' 2006). The heavy element mixture usually 

adopted is that of the Sun for population I stars, while for popu- 
lation II stars a given amount of a-element enhancement (for the 
same mixture) is taken into account. The precise values of the 
solar chemical abundances have become controversial because 
of the consider able impact of the new 3D hydrodynami c atmo- 
spheric models dAsplund et al.ll2005L l2009HCaffau et al.ll2010.) . 
We computed the stellar tra cks currently a v ailable in the 
database with the mixture of lAsplund et al.l (l2005h (here- 
after AS05) before the revised version of the same group 
(lAsplund et al.ll2009l hereafter AS09) was released; thus we as- 
sumed a reference value of solar photospheric metallicity of 
(Z/X)pj, - 0.0165. Since this is one of the main novelty of 
the present calculations with respect to other classical pre-MS 
models, we analysed the effect of the adopted solar mixture by 
comparing models computed adopting three different d istribu- 
tions, namely, the widely used iGrevesse & Noelsl (1 19931) (here- 
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Figure 5. Comparisons between non-grey models computed ad- 
opting the GN93, AS05 and AS09 heavy elements mixtures. The 
upper panel shows the mass range 1 .5-6.0 M©, while the bottom 
one the range 0.2 - 1 .0 M©. 



after GN93) and the two more recent AS05 and AS09. In each 
case, the same distribution was coherently used in both the opa- 
city tables (both the OPAL and the F05) and the computation of 
the nuclear-burning energy release. 

Figure |5] shows non-grey tracks computed with the quoted 
mixtures. Unfortunately, the atmospheric models adopted to ob- 
tain the boundary conditions are currently available for only one 
solar mixture (GS98 in CK03 and GN93 in BH05), hence the 
mixture is kept fixed in the atmosphere. 

As clearly shown in Fig. |5] the models computed with the 
two releases of the solar mixture of Asplund are almost indis- 
tinguishable. This is unsurprising because AS05 and AS09 are 
quite similar, with in particular almost identical global abund- 
ances in mass of both the iron group and a elements at a given 
metallicity Z. 

In contrast, the tracks computed with the GN93 mixture be- 
come dissimilar to the others when a sizeable radiative core 
develops and the stars leave the Hayashi line moving towards 
higher effective temperatures. In the ZAMS, the GN93 mod- 
els are hotter and brighter than the AS05 ones, any differ- 
ences becoming smaller and smaller as the mass decreases: 
ATeff ~ 300 K and AlogL/L^ ^ 0.02 for 6.0 Mq, ATeif * 90 K 
and A log L/Lq ^ 0.05 for 1 Mq, and negligible for 0.2 Mq. 

We note that the almost coincidence of the GN93 models 
with the others along the Hayashi track is the consequence of 
adopting the same heavy element mixture in the outermost lay- 
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Figure 6. As in Fig.|5]but for grey models. 



ers, whose radiative opacity governs the effective temperature of 
the Hayashi track, because the outer boundary conditions were 
computed with the same atmosphere models (see discussion in 
Sect.lza. 

Thus, to actually use different mixtures also in the atmo- 
sphere, we computed three additional sets of grey models by ad- 
opting the KS66 T(r) relationship with Tph = 2/3, even though 
grey boundary conditions should not be used in very low-mass 
stellar models (see e.g. Sect. I2.2l and references therein). 

Figure |6] shows the same kind of comparison as Fig. |5]but 
for the grey models. As one can see, the GN93 models in this 
case differ from the others also along the Hayashi track, as a 
consequence of having a different radiative opacity also in the 
outermost stellar layers that affects the outer boundary condition. 

For masses greater than 0.4 M© the GN93 models are the 
hottest along the whole pre-MS evolution. For these masses, an 
increase in the iron abundance causes a larger radiative opacity, 
thus a lower effective temperature and luminosity, exactly what 
occurs passing from GN93 to AS05 for the same metallicity Z. 
For less massive stars, the effect of updating the heavy element 
mixture from GN93 to AS05 on the pre-MS track is the result 
of two opposite trends: an increase in the H" contribution to the 
radiative opacity, as a consequence of the higher abundance of 
electron donors, and a decrease in the molecular opacity caused 
by the reduced oxygen abundance. 

The previous discussion refers to solar metallicity models, 
although we also performed comparisons for a lower metallicity, 
i.e. Z = 0.0002, for the three solar mixtures discussed above in 
the case of non-grey models. For metal-poor stars, we found that 
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the pre-MS tracks are insensitive to the heavy-element mixture 
in the HR diagram. 

2.5. Convection 

The boundaries of the convectively unstable regions are fixed 
by the classic Schwartzschild criterion; within these boundaries, 
we assume that the mixing timescale is short compared to the 
nuclear burning one. 

As is well known, one major shortcoming of stellar model- 
ling is the oversimplified treatment of convection, particularly 
in giant structures, such as RG and pre-MS stars, where large 
superadiabatic regions are present. The practical consequence is 
that we are not yet able to firmly predict the Teff and the radius 
of stars with an outer convective envelope, such as pre-MS ob- 
jects during the Hayashi track. The common app roach to stellar 
comp utation is to adopt the mixing length theory dBohm-Vitensd 
Il958l) for which the average convective efliciency depends on the 
mixing length £ - ffHp, where Hp is the pressure scale height 
and a is a free parameter that is calibrated with observations. We 
implemented th e mixing length scheme in the code following the 
prescriptions of ICox & Giuli ( 1968). 

The usual procedure for calibrating the mixing-length effi- 
ciency consists of constraining the a value by fitting the solar 
radius. For this reason, it might be useful to provide a brief de- 
scription of the main characteristics that a standard solar model 
(SSM) must satisfy. A SSM is a stellar model of 1 Mq such 
that at the age of the Sun (i.e. «; 4.57 Gyr) reproduces the 
observed solar luminosity, radius, and photospheric (Z/X)pi,. 
For this last quantity, we adopted the value 0.0165 from the 
heavy element mixture of Aspl und et alj (2005), the same mix- 
ture used in our reference models. In our SSM, the accuracy of 
the fit is very high, namely AL/L < 10"^, AR/R < 10 "*, and 
A(Z/X)pij/(Z/X)pj, < 4 ■ 10""*. The parameters to tune to achieve 
such a good fit are the initial metallicity Z, helium abundance 
Y, and the mixing-length parameter a. We note that the observed 
value of the photospheric (Z/X)ph does not represent the original 
one, since microscopic diffusion reduces the helium and heavy- 
element surface abundances, increasing the hydrogen one. 

Figure [T] shows the evolutionary tracks corresponding to 
three different SSMs provided by the current version of the 
FRANEC code (Yq = 0.2533, Zq = 0.01377). The physical in- 
gredients adopted to compute these models are exactly the same 
with the exception of the outer boundary conditions. In one case 
(dotted line), a simple grey atmosphere with the KS66 T(t) re- 
lationship is used, while in t he other two cases a non-g rey de- 
tailed atmosphere models by 'Brott & Hauschilda (|2005|) (solid 
line) and Castelli & Kurucz (2003) (dashed line) are adopted. 

To fit the present Sun, the three models require three different 
values of the MLT parameter, namely a - 1 .97 with the KS66 
grey T(t) relation, a = 1.68 and a = 1.75 with, respectively, 
BH05 and CK03 atmospheric modelfl As a result, the three 
models have a different pre-MS location and shape in the HR 
diagram, mainly near the Hayashi track, where the CK03 model 
is hotter than the BH05 one by about 60 K. The grey model is 
almost coincident with the CK03 one during the first contrac- 
tion, while it moves towards the BH05 track as the star evolves. 
Sim ilar results have akeady been discussed in the literatu re (see 
e.g.. lD'Antona & Mazzite"ililll994tlMontalban et al.ll2004l) . 
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Figure?. Comparisons between SSM evolutionary tracks com- 
puted adopting the different labelled boundary conditions, 
namely a grey atmo sphere plus the T(t) relationship by 
^Krishn a SwarnY| (|l966') (KS66, dotted line) and the non-grey at- 
mosphe res bvlBrot t & Hauschildt (2005) (BH05, solid line) and 
Castelli & Kuruca (2003) (CK03, dashed line) connected to the 
interior at Tph = 10. The symbol marks the position of the 
Sun. 



Although a better approach is still not possible, we recall 
that one should be careful when extending the solar-calibrated 
a to other stars because, as discussed in s everal papers, 
this could introduce systematic errors (see e.g., 'Pedersen et ajj 
1990; Ca nuto & Mazzitehi 1992; D'Antona & MazzitelH 1994 
Salaris & Cassisil 119961; iMo ntalban et al. 2004). We also note 
that there is no reason throughout the whole ev olution of the 
star to assume a constant value of a (see e.g., ISiess & Livid 
1 1 9971 ; [Piau et al.l |20 1 11 and references therein). In this respect, 
the progresses in 2D and 3D radiative hy drodynamical simu- 
lations of convection in the past decade ( Ludwig et alJ 119991 ; 



iTrampedach et al.lll999l ; iTrampedachI 12007 *) provided evidence 
of both variations in a, which can depends on both Teff and log g, 
and the possibility o f theoretically calibratin g the mixing length 
parameter (see also IMontalban & D'Antona .2006.) . Moreover, 
there appears to be a lower efficiency of the superadiabatic con- 
vection in low-mass stars (a ^ I) with respect to the solar- 
calibrated one according to observations of some binary systerns 
(see e.g., ISimon et al] 120001; ISteffen et al.i 2001; Stassun et al 
20041). and ^Li data of young o pen clusters (.Ventura et aL.199 



D'Antona & Montalbanll2003h . 



' Note that the mixing length parameter in the atmosphere a.,,tm differs 
from that used in the interior, namely a^,^ = 2.0 and a^in, = 1.25 for 
BH05 and CK03, respectively. 



This uncertainty in the mixing length efficiency also affects 
both mass and age estimates of young stars, since these are often 
obtained by comparing the stellar position in the HR diagram 
with the theoretical tracks, whose effective temperature depends 
on a whenever a convective envelope is present. 

Figure [8] shows the pre-MS tracks of a 0.8 Mq model with 
a = 1.90 and of four other masses in the range 0.8 - 0.95 M© 
with the solar-calibrated mixing length parameter, i. e. a = 1 .68. 
The points mark the intersection between the first tracks and the 
other four. As expected, the greatest differences occur near the 
Hayashi track at the point 'a' where the same luminosity and 
temperature of the 0.8 M© a = 1.9 model is reproduced by a 
0.9 Mq models with a - 1.68. The different choice of a trans- 
lates into differences in effective temperature of 90 - 100 K, 
which coiTespond to uncertainties of about 0.1 - 0.15 M© in the 
first phase of pre-MS contraction. This is consequently reflected 
in an uncertainty in age, because stars of different masses evolve 
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Figure 8. Figure shows that the same position in the HR diagram 
can be reached by stars of different masses if the a parameter 
is changed. This effect is shown with respect to a 0.8 M© star 
model with a - 1.9 (solid line) by comparison with models with 
a = 1.68 and masses in the range 0.8 - 0.95 M© (dashed lines). 
The boundary conditions are the same as our reference set of 
models. 



on different timescales; the largest effect on age occurring at the 
point 'b' where the relative difference is about 12%, whereas at 
the points 'a', 'c' is of the order of 6 - 7%. 

Given this intrinsic weakness in the mixing-length solar cal- 
ibration and the great impact of the efficiency of super-adiabatic 
convection, we decided to provide pre-MS models with three 
values of a, namely 1 .2, 1 .68 (standard solar model obtained us- 
ing the AS05 mixture and BH05 boundary conditions), and 1.9. 

2.6. Nuclear network and deuterium burning 

The current version of FRANEC code follows the burning 
of 26 elements, in particular the burning of light elements 
(D,^He, ''Li, ^Li, ''Be, and "B) starting from the pre-MS phase, 
with initial abundances of ■'He, ^Li, ^Li, ''Be, and "B from 
iGeiss & Gloecklen (Il998h . The cross-sections of the nuclear re- 
actions relevant to the hydrogen (pp chain and CNO cycle) 
and light element burning are assumed to be those of the 
NACRE compilation (Angulo et alj|1999h . with the exception 
of that of ''*N( p,y)^^0, which is from the LUNA collabo ration 



M 



(llmbriani et alfcoOStlLUNA Collaboration et al.ll2006allbh . 

While the burning of Li, Be, and B does not contribute much 
to the stellar luminosity because of their low abundances, deu- 
terium plays a crucial role during the first evolution of the pre- 
MS stars. The energy released by D-burning temporarily decel- 
erates the star's contraction, which accelerates again once this 
element has been exhausted. The observable effect of D-burning 
is to increase the stellar luminosity for a given age in the first 
few millions years (1 - 10 Myr, depending on the stellar mass). 

Recent estimate s based on the WMAP results (see e.g. 
iBennett et al.l 7005) and the standard model of cosmological 
nucleosynthesis, found a numerical abundance of the primordial 
deuterium relative to hydrogen (D/H)p = 2.75^q jg ■ 10"^ (see e.g. 
ICvburt et al]|2004lrl that corresponds to a fractional abundance 
in mass in the range 3.8 ■ 10~^ ^ ^d.^ ^■^. ' 10 ~^- A similar D- 
abundance was found bv lPettini et al.l (l2008h and lSteigman et all 
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Figure 9. Effects of two different assumptions about the initial 
deuterium abundance, namely Xd = 2 ■ 10"^ (dashed line) and 
Xd = 4 ■ 10"^ (solid line). Upper panel: time evolution of the 
luminosity for tracks in the mass range 0.2-6 M©. Lower panel: 
HR diagram for isochrones of ages 1, 2, 5, 10, 20, and 50 Myr. 



** By (D/H) we mean the ratio of the numerical abundances of deu- 
terium to hydrogen. 



( I2007h . On the other hand, for the lo cal ISM the value is lower , 
as a consequence of stellar astration. i Geiss & Gloecklen (Il998h 
found (D/H)isM ~ 2.1 ■ 10"^ (Xd ^ 3 • 10"^), while more re- 
cent observations give ( p/H)isM = 1-88 ± 0.11 • 10~^ that is 
Xn~2.6-ir^(se e also' Vi dal-Madiar et al Jl 1 9981: iLinskvl 1 998t 
iLinskv et al.ll2006l: Steigma neral.ll2007t) . 

Thus, Xd =4-10"^ and Xd = 2 ■ 10"^ should be suitable 
for, respectively, population II and population I pre-MS stellar 
models. 

Figure|9]shows pre-MS models (upper panel) and isochrones 
(bottom panel) computed with Xd = 4 ■ 10"^ and Xd = 2 ■ 10"^. 
As expected, the duration of the deceleration caused by the D- 
burning gets shorter with decreasing initial deuterium abund- 
ance. The extent of this effect depends on the mass of the star, 
since the deuterium exhaustion occurs at progressively earlier 
stages as the stellar mass increases. In particular, for stars with 
M <: 3.0 M0, the differences in the evolutionary times are limited 
to very young ages, well below 1 Myr Thus, in this mass range, 
the two sets of isochrones converge for ages older than 1 Myr, 
in particular for a 0.2 M© model, the effect of the adopted initial 
deuterium abundance is not negligible until much older ages, 
namely 10 Myr We also note that we discuss is the amount of 
deuterium actually left by the previous protostellar evolution in 
Sect.lZTl 

Since one of the novelties of the current set of models is 
the ^'*N(p,7')'^0 cross-section provided by the LUNA collabora- 
tion, a detailed discussion of the changes caused by replacing the 
older but still widely used NACRE cross-section with the new 
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Figure 10. Comparison between models with Z = 0.015, com- 
puted adopting the ''*N(p,7)'^0 reaction rate from LUNA (solid 
line) and NACRE (dashed line). Upper panel: evolutionary 
tracks in the HR diagram in the mass range 1-6 Mq. Bottom 
panel: isochrones for 1, 20 and 50 Myr 



one might be of some interest. The '^N(p,7)'^0 is the slowest 
reaction of the CNO cycle and thus governs the efficiency of the 
central H-burning in stars more massive than 1 .2 - 1 .5 M© (the 
exact value depending on the chemical composition), in which 
the CNO cycle is dominant with respect to the pp chain. 

Figure [10] shows a set of models and isochrones with 
Z - 0.015 computed adopting the NACRE (dashed line) and 
the LUNA (solid line) '''N(p,7')'^0 cross-sections. As expec- 
ted, the differences between the two sets of models become 
not negligible only for stars more massive than 1.2 - 1.5 M© 
when they approach the ZAMS, that is, when the CNO cycle be- 
comes a significant source of energy. The ZAMS models com- 
puted with the LUNA cross-section are hotter than those com- 
puted with the NACRE one, the former being smaller than the 
latter. The greatest differences in Tefr occur for higher masses 
(ATeff ~ 350 K for 6 Mq), while they decrease to 250 K for 
3 Mo and to 40 K for 1.5 Mq. 

Figure [TT] shows the ratio of the CNO to pp luminosities 
(LcNO and Lpp) versus the central temperature (Tc) for ZAMS 
models of different masses and metallicities computed adopt- 
ing the two quoted cross-sections. For a given value of the cent- 
ral temperature, the Lcno/^pp ratio drastically decreases as the 
metallicity becomes increasingly lower as a consequence of the 
lower C, N, and O abundances. Thus, the critical mass above 
which the effect of the adopted ''*N(p,7)'^0 cross-section on 
the ZAMS location is not negligible increases with decreasing 
metallicity. As an example, this critical mass rises from about 
1.5 Mo for Z = 0.015 to ^ 2 Mq for Z = 0.0002. We note that 
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Figure 11. Ratio of the CNO to pp luminosity (Lcato and Lp,,) as 
a function of the central temperature T^ for ZAMS models of dif- 
ferent masses, computed adopting the ''*N(p,y) '^O from NACRE 
(red lines) and LUNA (blue lines) for the labelled chemical com- 
positions. 



in the plane (log T^, log LcNo/^pp) the distance between the red 
and blue curves, computed respectively with NACRE and LUNA 
cross-section, is roughly independent of the metallicity. 

Regarding the ZAMS position in the HR diagram. Fig. [TO] 
the temperature differences decrease with metallicity: indeed for 
Z = 0.0002, we obtained ATes ~ 30 K for models of 2 M^, 
200 K for 3.0 Mq, and 300 K for 6 Mo. 

2. 7. Initial conditions 

The choice of initial conditions, i.e. the starting model, is im- 
portant when modelling the early phases of pre-MS evolution. 
The correct and physically consistent approach would be to take 
as a starting model the structure left at the end of the previ- 
ous hydrodynamical evolution of the protostar, when the main 
mass- accretion process is finished. In this respect, as early as 
1983 Stabler introduced the very useful concept of "birthline", 
defined as "the locus in the HR diagram where pre-MS stars of 
various masses should first appear as visible objects'" (Stahleg 
Il983l) . In agreement with this definition, realistic pre-MS mod- 
els should start from the birthli ne, which should also play the 
role of a zero-age isochrone (Pal la & Stahler|[l999l) . 

However, the hydrodynamical protostellar evolution is still 
largely debated and not yet settled. Thus, we decided to fol- 
low the standard app roach outlined by the pioneering contri- 
bution of'Iben' ('1965') and followed until now by many authors 
(Baraff^e et al. 1998; Siess et al. 2000; Yi et al. 2001; Potte r et all 
i2007l: iDi Criscienzo et al.ll20 09). which consists of choosing as 
an initial model a fully convective structure of very large radius 
and luminosity, i.e. a point along the Hayashi track, from which 
the star begins its quasi-static contraction at constant mass, i.e. 
neglecting any mass accretion episode. The problem consists 
of understanding whether this simplistic assumption signific- 
antly affects the whole evolution of pre-MS models or only the 
very early stages. Many works agree in predicting that after the 
end of the main mass-accretion phase, the evolution quickly 
converge to that of standard hydrostatic models ( Stabler 1983|; 
Palla & Stabler 199ll ll992l) . Ho wever, there is not yet complet e 
agreement with this conclusion dWuchterl & Tscharnuteil2003h . 
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Figure 12. Evolutionary tracks in the mass range 0.2-6.0 M© 
(dotted line) with superimposed 0.5 Myr (dashed line) and 1 Myr 
(so/ic/ line) isochrones. The birthhne by IStahler & Pallal (|2004|) 
is also shown (star-dashed line). 



We focus mainly on ages older than 1 Myr since the cur- 
rent theoretical predictions are generally considered unreli- 
able for younger ages (see e.g., iTout et al1ll999l : ISiessI l2001t 
LBaraffe et al. 2002 ). Figure [12] shows evolutionary tracks (0.2 - 
6.0 Mq) in the HR dia gram with the superimposed birthline of 
IStahler & Pallal (l2004 . and our isochrones of 0.5 and 1 Myr. 
Standard stellar pre-MS models intersect the birthline at different 
ages, depending on the mass, but almost always before 1 Myr 

Another point that deserves to be discussed is the deuterium 
abund ance of stars on the birthline. According to lStahler & Pallal 
(l2004l) . protostars start to burn deuterium once a total mass 
of 0.3 Mq is reached, if a constant mass accretion rate of 
10"^ Mq yr"' is adopted. As a consequence, the initial deuterium 
abundance available on the birthline at the beginning of the pre- 
MS phase is exactly the value present in the ISM for stars of 
mass lower than 0.3 Mq, while it decreases significantly with in- 
creasing total mass. According to Stabler & Palla (2004), stars 
with a mass higher than about 1 Mq should almost completely 
deplete the deuterium during their previous protostellar phase. 

In contrast, in standard pre-MS models, the initial deuterium 
abundance is usually the same as that of the ISM. In principle, 
such a difference could affect the early pre-MS evolution, but 
in practice the effect is negligible, since in standard models the 
deuterium burning onset occurs at luminosities much higher than 
those corresponding to the birthline and when the model reaches 
it the deuterium depletion agrees with that provided by protostel- 
lar computations. As a test, we compared the evolution of a 
0. 9 Mq model provided by our standard computation with that 
of IStahler & Pallal (l2004l) . which takes into account protostellar 
accretion. Starting from the same interstellar deuterium abund- 
ance, our model attains the luminosity of the birthline experien- 
cing the same amount of deuterium depletion, namely a tenth of 
the interstellar value, as the model emerging from a protostellar 
evolution. 

The definition of proper starting model is also of interest be- 
cause the age issue. Usually, the age assigned to pre-MS stars by 
means of standard models is simply the interval of time required 
for the model to attain the observed radius and luminosity start- 
ing from the arbitrary initial conditions of high luminosity and 
large radius. This inferred age is inaccurate, since it completely 
neglects the protostellar phase, an inaccuracy that grows as the 
age of the star decreases. The correct approach should be to add 



the protostellar ag e to the pre-MS one, measured from the birth- 
line (IStahleill983l) . On the other hand, this inaccuracy should be 
between O .I Myr and 0.6 M yr, the time needed to form a star ac- 
cording to lStahler & Pallal (l2004l) models which adopt a constant 
mass accretion rate of I0~^ M© yr"'. 



3. Comparison with different authors 

We compare our reference set of tracks with some of the 
most recent and used pre-MS m odels avail able in the hterat- 
urefl na mely.lBaraffeetal.l (fT99| BCAH98), 'Siess et al." (I20Q3 
SD OOI. lYi et all (l200lL YYOl ), Potter et al. (2008, DSEP08), 
and lDi Criscienzo et al.! ("200?, DVD09E1 

The upper panels of Fig. [13] [14] and [15] show, respectively, 
the evolutionary tracks of 0.4 Mq, 1.0 Mq, and 3.0 Mq with 
Z «; 0.02. Table [2] lists the main characteristics (EOS, 
radiative opacity, boundary conditions, convection scheme, 
cross-sections, and initial deuterium abundance) of the quoted 
models. 



M = 0.4 Mq 

Figure [13] (upper left panel) shows the HR diagram with the 
evolutionary tracks of 0.4 Mq. The models by BCAH98 with 
a = 1, DVD09 and FRANEC have a quite similar location in 
the HR diagram for log L/Lq < - 1 as a result of using similar 
EOS and boundary conditions. Figure [13] (bottom left panel) 
shows that the effect of adopting different a values becomes 
progressively negligible as the star approaches the ZAMS. 
More quantitatively, changing a from 1.68 to 1.2 produces a 
difference in effective temperature of about 150 - 170 K for ages 
of 1 Myr, difference that decreases to 30 K at logL/Lg = -1 
and to 10- 15 Kin ZAMS. 

The model of DSEP08 is slighdy colder (about 40 K) than 
those of BCAH98, DVD09, and FRANEC as it approaches the 
ZAMS; this difference is not easy to understand because the rel- 
evant input physics are similar Moreover, we checked that the 
difference of a between DSEP08 and our models would produce 
a ZAMS hotter by about only 3 K (see bottom left panel of Fig. 
[T3l l. The choice of Tph could account, at least in part, for the 
observed discrepancy. DSEP08 uses Tph such that T(rph) - Teff, 
which leads to a variable value often lower than one. We verified 
that changing Tph from 10 to the value such that T(Tph) = Teif pro- 
duces ZAMS about 15 - 20 K colder We also note that DSEP08 
model has a lower original metal and helium abundances, but the 
two effects cancel each other. 

The tracks of SDOO and YYOl show the greatest differences 
in both morphology and position on HR diagram compared to 
th e others. SDOO u se a modified version of the EOS described 
in Pols e t all d 1 995b that, as discussed in Sect. 12.11 should pro- 
duce tracks colder than ours of about 35 K. A further shift (about 
5 - 10 K) towards lower effective temperatures due to the lower 
helium abundance is expected. However, Fig. [13] shows that the 
track is hotter than the other models for log L/Lq < -0.5, reach- 
ing in ZAMS ATeff ~ 150 K. This result is not easy to explain. 
One might speculate that a not negligible role in the deviation 



' The models selected for the comparisons obviously do not 
exhaust t he sample of pre-MS tracks a vailable in the l it erature , 
see e.g. 'D' Antona & Mazzitellil ( 119971) , ICharbonnel et alj J1999I) , 
Palla & Stabler ( 1 99^^ 

'" DVD09, which is the updated version of lD'Antona & MazziteUj 
( 1997), deal with low metallicities; the models with Z = 0.02 were com- 
puted by Marcella Di Criscienzo with the same evolutionary code de- 
scribed in DVD09. 
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Table 2. Summary of the main characteristics of the selected grid of models. The columns provide the adopted EOS, radiative 
opacity (the heavy element mixture is specified), boundary condition characteristics, superadiabatic convection treatment, nuclear 
cross-sections and initial deuterium abundance. All these codes use a solar-calibrated convection efficiency, except DVD09. 



Code: 


EOS 


Radiative opacity 
(Mixture) 


Boundary conditions 


Convection 


Reaction rates 


Xd 


FRANEC 


OPAL06 


OPAL; 

F05 
(AS05) 


non-grey, Tph = 10 
Brott & Hauschildt (2005^ 


MLT 
a= 1.68 


NACRE; 

LUNA 

(i'*N(p,r)i5o) 


2 ■ 10"\ Z > 0.008 
4 ■ 10-\ Z < 0.008 




ffatm = 2; 

Castelli & Kurucz (2003) 

ffatm = 1-25 




BCAH98 


SCVH95 


OPAL; 

AF94 

(GN93) 


non-grey, Tpi, = 100 
Allard & Hauschildt (1997) 

C*'atm = 1 


MLT 
0-= 1.9" 


CF88 


2 ■ 10"' 


DSEP08 


Chabover & Kim ri995'); 


OPAL; 

F05 
(GS98) 


non-grey, Trp^=Teff 


MLT 
a = 1.938 


AD98; 

LUNA04 
(l*N(p,7)'50) 







Q^atm — 2; 

Castelh & Kurucz (2003) 

ffatm = 1.25 




DVD09 


OPAL05; 
SCVH95 


OPAL; 

F05 
(GS99) 


non-grey '' , M < 2 M© 
Allard & Hauschildt (1997). t„,, = 3 

Q'atm = 1; 

Heiter et al. (2002). FST. T„h = 10; 

grey, M > 2 Mq 


MLT'- 

a = 2.0; 

FST 


NACRE 


4 ■ 10^5 


SDOO 


modified version of 
PTEH 


OPAL; 
AF94 
(GN93) 


non-grey, Tph = 10 
Kurucz (1991); Plez (1992) 


MLT 
a= 1.6 


CF88 


2 ■ 10-5 


YYOl 


OPAL96; 
Chabover & Kim (1995^ 


OPAL; 

AF94 
(GN93) 


grey atmosphere, 

Tph = 2/3 


MLT 

a = 1.743 


BP92 


(...) 



No tes. '°' BCAH98 m odels have been computed also for a = I; we used this latter val ue of a for the comparisons of 0.4 M o- DVD09 adopt 
the lAllardetalJ ( 119971) atmospheric models for M< O.6M0 and the lHeiteret alj J2002I) ones for M> O.6M0. '•"'> DVD09 models were computed 
adopting the MLT formalism for M< 0.6 Mo and the FST for M> 0.6 Mq. 



Refer ences. PTEH. iPols et all (119951) : SCVH95,ISaumon et al.l (11995'); AD9 8. lAdelbereer et al J d 1998h : B P92.lBahcaU & Pinsomieaultl (|1992h : 
CF88. ICaughlan & Fowled dl988b: LUNA04.l!m briani et al. (2004); GS99, Gre vesse & Sauvall ( ll999l) ;"AF94. lAiexander & FergusoiJ ( 119941) ; FST. 
Full Spectrum of Turbulence dCanuto & Mazzitelliil991i ; iCanuto et al.iil996) . 



of the SDOO track is played by t he scheme adopted to integrate 
the atmosphere. As discussed in ISiess et al.l ( l2000l) . the bound- 
ary conditions are obtained by integrating the T(t) relationship 
resulting from a fit of atmospheric models. 

To test such a working hypothesis, we computed additional 
tracks (shown in the bottom right panel of Fig. [T3] ). using two 
different T(t) relationships, namely the classic KS66 (dashed 
line) and that obtained by interpolating the BH05 atmospheric 
models (dotted line), the latter representing the closest approx- 
imation to the SDOO choice. Both the models were computed 
using the same low-temperature radiative opacity adopted by 
SDOO, i.e. AF94. The dotted line model is hotter than the ref- 
erence one (solid line) of about 50 - 60 K in ZAMS. Moreover, 
as discussed in Sect. 12.31 the integration of a hydrostatic atmo- 
sphere for 0.4 M0 stars is very sensitive to the extrapolation 
method used for the low-temperature opacity tables. This effect 
can contribute to increasing the differences in effective temper- 
ature. 

For YYOl, the track morphology is completely different 
from the others; the ZAMS is 500 - 550 K colder than the 
BCAH98, DSEP08, DVD09, and FRANEC ones. This large dis- 
crepancy can be justified by neither the different initial helium 
abundance nor the adoption of the AF94 low-temperature radiat- 
ive opacity in concomitance of the use of a grey T(t) relationship 



(see bottom right panel of Fig. [T3] ). It is not easy to unambigu- 
ously identify the reason for this behaviour without a full evolu- 
tionary computation with the same input physics. However. lSiessI 
(1200 ih did prove that the location in the HR diagram of pre-MS 
tracks of low-mass stars is very sensitive to the contribution of 
molecular hydrogen in the EOS, showing that models computed 
with an EOS that does not account for H2 are significantly less 
luminous (AlogL/Lg ^ 0.6 for a 0.3 M©) and colder (about 
1000 K for a O.3M0) than those computed with an EOS that 
does take it into account. This seems to be the case for YYOl, 
because for temperatures lower than 10'' K, in the regions not 
covered by the OPAL EOS96, they adopt the Saha equation that 
accounts for a single state of ionization of atomic hydrogen and 
metals plus double states of ionization of he lium but they neglect 
the contribution of H2 (see lGuenther et aflri992i) . 

The right upper panel of Fig. [T3] shows the temporal evolu- 
tion of the luminosity. The models by BCAH98, DSEP08, and 
FRANEC agree for logt(yr) > 6.5 - 7. The non-negligible dis- 
crepancy at earlier stages between the FRANEC and DSEP08 



models (AlogL/Lg 



0.3 - 0.4 at logt(yr) = 6) is caused 



by the different treatments of D-burning, DSEP08 assuming a 
zero initial deuterium abundance. 

For log t(yr) » 6, BCAH98 is significantly fainter than 
SDOO and FRANEC, although these tracks adopt the same ini- 
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Figure 13. Upper panels: comparisons of 0.4 M© tracks from different authors listed in Table|2]for Z ^ 0.02, in the HR diagram and 
in the (logt(yr), logL/L©) plane. Bottom left panel: 0.4 M© models computed with three values of the mixing length parameter, 
namely a - 1 .2, 1 .68, and 1 .9. Bottom right panel: 0.4 M© models with three different choices of the adopted boundary conditions: 
non-grey models as described in Sect. \2.2\(solid line), KS66T(t) relation plus AF94 low-temperature opacity tables (dashed line) 
and T(t) relation interpolated from the iBrott & Hauschilda (l2005h atmosphere models plus AF94 low-temperature opacity tables 
(dotted line). 



tial deuterium abundance. This behaviour can be related, at 
least in part, to the fact that BCAH98 track starts its evol- 
ution from an initial model with a central temperature high 
enough to ignite deuterium, at variance with SDOO and FRANEC, 
which start their evolution at m uch lower central temperature. 
iD'Antona & MontalbanI (l2006l) showed that the initial central 
temperature significantly afi'ects the D-burning phase, in particu- 
lar tracks evolved from hotter initial models are fainter at a fixed 
age, as in the case of BCAH98. However, this effect is limited to 
a few Myr for a 0.4 Mq, which is the typical deuterium burning 
timescale. 

The DVD09 model displays a quite peculiar behaviour that 
cannot be simply explained by the different initial deuterium 
abundance (Xd = 4 ■ 10"^). For logt(yr) < 7.5, it is sys- 
tematically less luminous than the other models, with a max- 
imum discrepancy in luminosity of about AlogL/Lg ^ 0.9 at 
1 Myr One might speculate that such a large difference in lu- 
minosity for such a low-mass star is related to the D-burning 
phase at the beginning of the Hayashi track that is n ot in- 
cluded in DVD09 computations. iDi Criscienzo et al.l (1200 9) dis- 
cussed how, owing to the lack of atmospheric structures for 
logg(cms"^) < 3.5, they had to skip the deuterium burning 
during the early evolutionary phases. The DVD09 0.4 M© track 



begins at logg(cms"^) w 4.00, at a luminosity lower than the 
other models (logL/Lg « -1), BCAH98 included, who use the 
same atmospheric models. As the star approaches the ZAMS the 
DVD09 model eventually converges to BCAH98. 

As already seen in the HR diagram, the YYOl model shows 
a temporal evolution of luminosity quite different from the other 
ones. While the other models converge after the D-burning 
phase, the YYOl one becomes progressively less luminous 
until it reaches the ZAMS with a maximum difference of about 
AlogL/Lg^ 0.2. 

M = 1.0 Mo 

Figure [14] compares the tracks of 1 Mq. All the models predict 
almost the same luminosity for the ZAMS. The effective 
temperatures of the models of BCAH98, DVD09, YYOl, and 
FRANEC are still in reciprocal good agreement, whereas DSEP08 
and SDOO ZAMS are, respectively, hotter by about 94 K and 
colder by about 210 K than ours. 

As previously shown in Sect. 12.11 the ZAMS position of 1 
Mq stars in the HR diagram is quite insensitive to the EOS, 
i.e. the difference between PTEH and OPAL models is about 
15 K. In contrast, as shown in Sect. 12.31 the chosen heavy- 
element mixture has a significant effect on the ZAMS, in par- 
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Figure 14. Upper panels: comparisons of 1 M© tracks from different authors listed in Table |2] with Z ^ 0.02, in the HR diagram 
and in the (logt(yr), logL/L©) plane. Bottom left panel: 1 M© models computed with three values of a, namely a - 1.6, 1.68, 
and 1.9 for Z - 0.02 and Y - 0.288. Bottom right panel: 1 M© tracks with three different original helium abundances, Y = 0.270, 
Y = 0.280, and Y = 0.288 with Z = 0.02. 



ticular substituting the old GN93 with the most recent AS05 
mixture leads to colder (a; 90 K) and fainter (a; 0.05) 1 M© mod- 
els. DSEP08 adopt GS98 and therefore the difference should be 
smaller However, for DSEP08 this shift in the effective temper- 
ature is largely balanced by the different boundary conditions 
chosen. We verified that when we adopt the same prescription 
as DSEP08 (T(Tph) = Tes) the position of the ZAMS is shifted 
towards lower temperatures by about 50 K with respect to the 
model computed with Tph = 10. Moreover, the DSEP08 model 
adopts ff a 1.9, which would produce a ZAMS hotter than that 
computed with a - 1.68 for about 70 K, as shown in the bot- 
tom left panel of Fig. [14] Therefore, the net effect of the differ- 
ent heavy-element mixture, boundary conditions, and a values 
between the FRANEC and DSEP08 models is to produce a differ- 
ence in the effective temperature on the ZAMS of the order of 
100 K, which is fully compatible with the 94 K actually present. 
The slight difference between the initial metal and helium abund- 
ances of DSEP08 and our model is inconsequential because the 
effect of increasing helium is counterbalanced by the concomit- 
ant reduction in metals. 

The track of SDOO is computed with a - 1.6, which should 
produce a ZAMS colder than ours by about 40 K (bottom left 
panel of Fig. fT4l) . In this case, the total effect of the different 
heavy-element mixture and mixing length parameter would be a 
ZAMS model hotter than the FRANEC one of about 60 K, but ac- 



tually colder of 2 1 K. This large discrepancy can not be entirely 
justified by the different helium abundance. As shown in the bot- 
tom right panel of Fig. [141 even a difference larger than that 
between the FRANEC and SDOO models, namely AY = -0.018, 
would reduce the ZAMS effective temperature by about 100 K 
making SDOO colder than our track of about 40 K. However, the 
peculiar behaviour of th e SDOO 1 Mq model wi th Z - 0.02 was 
previously discussed bv lMontalban et aP (l2004l) . 

The maximum difference in temperature along the Hayashi 
track reaches almost 300 K between DVD09, which is the hottest 
model, and SDOO, the coolest. The agreement among the other 
tracks is better, i.e. ATeff ~ 60 K between FRANEC and BCAH98. 
With the exceptions of DVD09 and YYOl, the pre-MS models 
tend to converge at the end of the Hayashi track when stars move 
towards the ZAMS. We note that the lower the helium abundance 
the colder the Hayashi track, with a shift in effective temperature 
of about 36 K for AY = 0.018 at logL/L^ = 0.5 (see Fig.[ll 
bottom right panel). However, the dominant role in determining 
the position of the Hayashi track is played by the efficiency of the 
superadiabatic convection. Figure[T4lbottom left panel shows the 
effect of changing the a value between 1 .9 and 1 .6 in the interior 
(ATeff ~ 150 K at logL/L^ = 0.5). 

The upper right panel in Fig. [14] shows the time evolution 
of luminosity. At the age of 1 Myr, the maximum difference 
between the models is AlogL/Lg ^ 0.2 and decreases as the 
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Figure 15. Upper panels: comparisons of 3 M© tracks from different authors listed in Table|2]for Z ^ 0.02, in the HR diagram and in 
the (log t (yr), log L/L©) plane. Bottom left panel: 3 M© models with three initial helium abundances, namely Y = 0.270, Y - 0.280, 
and Y - 0.288 with Z = 0.02. Bottom right panel: 3 Mq tracks computed adopting the LUNA '^N(p,y)'^0 reaction rate {solid line) 
and the NACRE one {dashed line). 



models evolve, reaching A log L/L^ k 0.03 - 0.04 in ZAMS. 

M = 3.0 Mo 

Figure [15] shows the evolution of a 3 Mq. The models by 
BCAH98 are unavailable for M > 1.4 Mq, while the track 
of YYOl has only a few points making the comparison quite 
difficult and not so meaningful. Thus, we present only the 
comparison with DSEP08, DVD09, and SDOO ti-acks. Moreover, 
we do not discuss the location of the Hayashi track, since 
the star takes about 1 Myr to leave it, but focus instead on 
describing the discrepancies between the models for two 
specific evolutionary phases, namely the first relative maximum 
in effective temperature before the steep drop in luminosity 
(hereafter hoolrT) and the ZAMS. We note that these two points 
are insensitive to the efficiency of the superadiabatic convection 
owing to the lack of a significant convective envelope. 

The discrepancy in the hook effective temperature between 
the models is about 310 K for DSEP08 and FRANEC, 360 K for 
DVD09 and FRANEC, and about -150 K for SDOO and FRANEC. 
To help us understand the origin of these differences, the bot- 
tom panels of Fig. [15] show models computed with three values 
of Y namely 0.288, 0.280, and 0.270 for Z = 0.02 (bottom left 



" This phase corresponds approximately to the first model which is 
completely supported by central hydrogen burning with the secondary 
elements not yet being in equilibrium. 



panel), and with two values of the '^N(p,7)'^0 nuclear cross- 
section (bottom right panel), at fixed chemical composition, one 
by the NACRE compilation and the other by the latest release of 
the LUNA facility. The original helium abundance variation ac- 
counts for differences with respect to the FRANEC model of about 
-150 K for DSEP08 and -100 K for DVD09, whereas the effect 
of the '"^N(p,7)'^0 reaction rate becomes important only near 
the ZAMS and is completely negligible in the hook. We note 
that DSEP08 is also slightly metal poorer than the other models: 
this leads to a shift of about 1 10 K towards higher effective tem- 
peratures, hence the effect of the lower initial helium abundance 
is approximatively counterbalanced by the lower metal content. 
We recall that the position of the hook is also quite sensitive to 
the heavy-element mixture, that, as discussed in Sect. 12.41 would 
account for a difference of about 300 K in effective temperature 
and about 0.03 in log L/Lq between tracks computed with the 
GN93 and AS05 mixtures. A small correction to this estimate 
is expected, since DSEP08 and DVD09 use a slightly different 
heavy-element mixture. Thus, because of the effects of the ad- 
opted chemical composition and solar mixture, we would expect 
a net shift towards higher effective temperatures than our model 
of about 260 K for DSEP08 and 200 K for the DVD09, which 
is not enough to justify the observed differences of, respectively, 
310 K and 360 K. 
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Figure 16. Upper panels: comparisons between DSEP08, YYOl, and FRANEC for Z a 0.006 in the HR diagram and in the (log t(yr), 
log L/Lq) plane for 0.4 M©, 1 M©, and 3 Mq models. Bottom panels: as in the upper panels but for a lower metallicity, namely 
Z * 0.0002. 



At the ZAMS location, we found discrepancies in effective 
temperature with respect to FRANEC of about 200 K for DSEP08, 
100 K for DVD09, and -400 K for SDOO. The effect of the differ- 
ent helium abundances and heavy element mixtures are similar 
to those mentioned above. In addition, as discussed in Sect. 12.61 
and shown in the right bottom panel of Fig. [15] the ZAMS loc- 
ation is very sensitive to the adopted '"^N(p,y)'^0 reaction rate. 
The difference in effective temperature between tracks computed 
with the NACRE and LUNA cross-section is about -230 K, a 
value that should also be representative of the difference between 
LUNA and Caughlan & Fowler (1988) reaction rates used by 
SDOO (see e.g.. iAngulo et al.il99 9). Owing to the different initial 
helium abundances and heavy-element mixtures adopted, we ex- 
pect the DVD09 and SDOO models to be respectively 30 K and 
60 K colder than the FRANEC one, while DVD09 in contrast is 
100 K hotter and SDOO is about 340 K colder than our predic- 
tion. 

DSEP08 use th e old LUNA ''*N(p,r)'^0 reaction rate 
dlmbriani et al.ll2004l) . which is about 6% greater than the latest 
LUNA release adopted in FRANEC. This difference would make 
the DSEP08 model colder than ours by about 20 K. Adding 
to this the other differences in the physical inputs adopted by 
DSEP08 with respect to FRANEC (i.e. initial helium and metal 
abundances and mixture) would account for a discrepancy of 
about 240 K, which is approximatively 40 K higher than what 
is observed. 



The pre-MS time evolution of the luminosity for this mass 
is shown in the right upper panel of Fig. [15] Before the hook 
(6 < logt(yr) < 6.4), FRANEC and SDOO models, for a fixed 
age, are less luminous than DVD09 (AlogL/Lg ^ 0.07 and 
AlogL/L0 K 0.1 respectively). The ZAMS age of DVD09 and 
FRANEC are in reciprocal good agreement ( At/t % 5%), whereas 
SDOO predict a ZAMS age older by about 20-30% than FRANEC 
and DVD09. 

In addition the luminosity evolution as a function of age of 
DSEP08 is in good agreement with our model, but we note that 
there is a peculiar bump near the ZAMS, which none of the other 
models show. 



3.1. Comparison between metal poor models. 

Figure [16] shows the comparison for sub-solar metallicity, 
namely Z x 0.006 (top panels) and Z x 0.0002 (bottom pan- 
els), between our models and those of DVD09, DSEP08, and 
YYOl[]j. We note that although the initial metallicities are very 
similar, if not identical, the initial helium abundances are quite 
different. Thus, the observed discrepancies in the HR diagram 
location, in particular near the ZAMS, are not only the result of 



'^ BCAH98 and SDOO metal-poor models are unavailable. We also 
note that DVD09 low-metallicity models are available only for M< 1.5 
Ms. 
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different physical inputs but also of the adopted initial chemical 
compositions. However, the models are in reasonable agreement 
with the exception of the YYOl 0.4 M© track, which is signi- 
ficantly colder than the others (by more than 300 K in ZAMS), 
probably as a consequence of the adopted EOS, as previously 
discussed. We note that, at variance with the case for Z = 0.02, 
the effect of adopting different heavy-element mixtures is small 
for models with Z ^ 0.006 and completely negligible in the case 
of Z «; 0.0002, as we discussed at the end of Sect. 12.41 

The agreement in the HR diagram is very good between 
FRANEC and DSEP08 across the whole range of masses for both 
the metallicities, with the exception of 0.4 M© which is slightly 
colder than ours by about 40 - 50 K. A quite good agreement 
is achieved also between FRANEC and YYOl for Z = 0.0054 
(0.4 Mo excluded), whereas the Z = 0.00017 YYOl tracks have 
a ZAMS that is colder than ours by about 3 - 4%. When com- 
paring DVD09 and FRANEC, very small differences are present 
near the Hayashi track for both the metallicities, which increase 
near the ZAMS of 1 .0 Mq where DVD09 models ai-e colder than 
DSEP08 and FRANEC by about 130 - 150 K. 

As one can see in the right panels of Fig. [16] the temporal 
evolution of the luminosity for metal-poor models displays the 
same features as discussed for the Z « 0.02 case. 
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Figure 17. HR diagr am of the core of the cluster IC 348 by 
iLuhman et al.l ( 1199 8') with our theoretical pre-MS track super- 
imposed for the mass range 0.2 Mq to 7.0 M© (dotted lines) and 
isochrones of ages 1 Myr (dashed line), 2 Myr (dot-dashed line), 
5 Myr (dot-dot-dashed line), and 40 Myr (solid line). 



4. Comparison with observations 

The above discussion confirms a well-known feature of pre-MS 
tracks and isochrones: their position in the HR diagram is very 
sensitive to many still poorly constrained factors. The corres- 
ponding uncertainty in the predicted Teff and L directly propag- 
ates into an uncertainty in the inferred masses and ages of ob- 
served young stars and consequently in the IMF and SFR of 
young stellar systems. Hence, it is of crucial importance to com- 
pare the theoretical predictions with observations. 

Figure [17] sho ws the HR dia gram of the core of the cluster 
IC 348 studied bv .Luhman et al.l (11998.) with our theoretical pre- 
MS tracks and isochrones superimposed. With the exception of 
the faintest tail, the agreement is quite good. 

However, a much more severe test of pre-MS models is bin- 
ary stars, especially detached, double-lined, eclipsing systems, 
which allow us to calibrate the theoretical tracks against stars 
of known mass. Figure [TS] shows th e comparison between the 
RX J0529.4+0041A binai-y system ( Covino et al.l 12000'. '2004^) 
and the present theoretical tracks a nd i sochrones. This binary has 
been e xtensively studied (see e . g.,lD' A ntona 2001; Baraffe et al. 
I2002t iMontalban et"an l2004t iMontalban & D' Antona. i2006) 
and provides a very useful and tight benchmark for pre-MS 
evolutionary tracks, because the two low-mass components 
(Mi= 1.27 + 0.01 Mo and M2= 0.93 + 0.01 Mq) are still in 
the pre-MS phase. Even if an exhaustive empirical test of the 
tracks should rely on a detailed discussion of the effects of the 
various sources of uncertainties, not least the chemical compos- 
ition, we can safely claim that the agreement is quite good and 
very encouraging. 

We note that among the current theories of binary formation 
there is not yet consensus on the coevality of the two stellar com- 
ponents, thus we should not in principle expect the two compon- 
e nts of a binary sy stem to have the same age (see e.g. comment 
in lPalla & Stahleil[200L) . This means that we should not pretend 
that both the stars of a binary are fitted by the same isochrone. 
On the other hand, the comparison between theory and observa- 
tion seems to suggest coevality, that is the compon ents of bin- 
aries form at about the same age (see e.g. Hartigan et al.lll994t 
iBrandner & Zinneckedll997blPalla & Stdiler.,2001.) . 
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Figure 18. HR diagram o f the binary system RX 
J0529.4-1-0041A dCovino et all |2004 with our theoretical 
pre-MS tracks of labelled mass (solid lines) and isochrones of 7 
Myr (dotted line) and 10 Myr (dashed line) superimposed. 



5. The database 

As is well known, the location in the HR diagram of pre-MS 
tracks, in particular the Hayashi line, strongly depends on the 
chemical composition, and particularly the metallicity. As an ex- 
ample of how large this effect can be. Fig. [19] shows 0.8 Mq 
models with different chemical compositions. Figure l20l shows 
the results for two selected metallicities, namely Z - 0.005 and 
0.008, corresponding to a difference in [Fe/H] of about 0.2 dex, 
roughly twice as much as the typical spectroscopic uncertainty. 
The shift in effective temperature of the Hayashi track at fixed 
mass is about 100 K. This means that to reproduce the Hayashi 
track location of a model with M= 0.8 Mq and Z - 0.008 with 
tracks of Z = 0.005, a decrease in the mass of 0.1 M© is required. 
Furthermore, this metallicity error also affects the age, since the 
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Figure 19. HR diagram with 0.8 Mq pre-MS tracks with the dif- 
ferent labelled chemical compositions. 
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Figure 20. Pre-MS tracks of 0.8 and 0.7 Mq with Z = 0.005 and 

Y = 0.258 (dashed lines) and of 0.8 Mq with Z = 0.008 and 

Y — 0.265 (solid line). 



contraction timescale of a 0.7 M© star is longer than that of a 
0.8 Mq starl This simple example gives an idea of the error in 
the inferred mass and age of stars in the pre-MS phase caused 
by an eiTor in the assumed chemical composition. Thus, when 
comparing data with theoretical pre-MS tracks, one should be 
very careful to use models with the same chemical composition 
(in particular the same metallicity) of the observed stars. 

This is why we provide a database of models with a very 
fine grid of Z and Y values. The present models are available 
for 19 metallicity values, from Z = 0.0002 to Z = 0.03. For each 
value of Z, we computed models for three different initial helium 
abundances Y, following the relation 



AY 



(1) 



where Yp represents the cosmological '*He abundance and 
AY/AZ is the Galactic helium-to-metal enrichment ratio. 



At logL/Lo=-0.4, where the model with Z = 0.008 starts to leave 



the Hayashi track, the age of the 0.7 Mg 
one of 0.8 Mn with Z = 0.008. 



model is 70% older than the 



For the cosmological value of Yp, we used both the re- 
cent WMAP estimate Y^ = 0.2485 (ICvburt et al.l l2004t 
Steigman 2006; Peim bert et al. 2007a,b) and an older estim- 
ate o f Yp = 0.230 dLequeux et all 119791: iPagel & SimonsonI 
119891; lOlive et all 1 19911) . usually referred to as canonical Yn 
in several stellar isochrone databases (e.g, [V^n denBerg et al] 
'2000'; 'Girardi et alj 120001: lYi et al.l l200ir iGirardi et aLlGOOlT 
Cariulo et al. 2004;). 

The estimated value of helium to metal enrichment is af- 
fected by several uncertainties, thus we decided to choose both 
the typical value AY/AZ = 2 and a hi gher value AY/AZ = 5 , 
whi ch seems to be the upper extreme (Pagel & Portinaril ll998l; 
Jimenez et al.ll2003HFivnn 2004: Gennaro et al. 2010). For each 
value of Yp, AY/AZ, and Z, we computed models in the mass 
range 0.2-7.0 M© with different mass steps: 0.05 M© in the mass 
range 0.2-1.0 Mq, 0.1 Mq in the range 1.0-2.0 Mq, 0.2 Mq for 
the mass range 2.0-4.0 Mq, and 0.5 Mq in the range 4.0-7.0 Mq. 

The models are available for three different values of the 
mixing length parameter a, namely 1.68 (solar calibrated), 1.2, 
and 1.9. 

All models adopt an initial deuterium abundance 
Xd - 4 ■ 10"'', which is a good estimation at least for 
population II stars. As discussed previously in Sect. 12.6! popu- 
lation I stars should have a lower deuterium abundances, hence 
models with Xq = 2 ■ 10"^ for Z > 0.008 are also available. 

The grid of models available in the database are summar- 
ized in Table |3] We also provide an additional set of models 
corresponding to the chemical composition of our SSM, namely 
Z = 0.01377 and Y = 0.2533. 

Isochrones for ages from 1 Myr to 100 Myr are available; the 
age spacing is 1 Myr for ages 1-20 Myr and 5 Myr between 
20 Myr and 100 Myr. 

Figures |2T] and |22] show a few examples of complete sets 
of evolutionary tracks and isochrones, respectively, as extracted 
from our database in the (logTeff, logL/Lg) theoretical plane. 
Figure |2T] also shows the isochrones of 1 Myr, which, as dis- 
cussed in Sect. 12.71 represent a first approximation of the birth- 
line; the dashed part of the tracks correspond to ages younger 
than 1 Myr, those more affected by theoretical uncertainties. 

Models and isochrones are already available in the theor- 
etical plane (logTefj, logL/Lg) and will be available in sev- 
eral photometric systems in the next future. The files of each 
track list: model number, log t (yr), central abundance by mass 
of hydrogen, logL/Lg, logTeff, central temperature and dens- 
ity (logTc and logpc, respectively), mass of the convective core 
Mcc, contribution to the total luminosity of the proton-proton 
chain (Lpp/Ltot) and CNO cycle (LcNof^tot) burning and of the 
gravitational energy (Lc,/Ltot), while the files of the isochrone list 
log L/Lq, log Teff , and the mass of the model (M/M©). 

6. Conclusions 

A growing amount of data of young stellar systems has prompted 
renewed interest in modelling the initial phase of stellar evolu- 
tion. 

The present paper describes a new set of pre-MS tracks and 
isochrones, which relies on state-of-the-art input physics (EOS, 
radiative and conductive opacity, atmospheric models and nuc- 
lear cross-section). To provide the astronomical community with 
a useful and versatile theoretical tool for the interpretation of 
observational data, the models have been computed for a very 
large and fine grid of chemical compositions. We have evaluated 
model data for 19 metallicities, ranging from Z - 0.0002 to 0.03, 
and three different helium abundances for each Z. For a fixed 
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Figure 21. Theoretical tracks in the mass range 0.2 - 7 M© for labelled chemical composition with a = 1 .68 and Xd - 4 ■ 10"^ with 
superimposed the corresponding isochrone of 1 Myr (dashed line). The dotted line represents the evolutionary phases with ages 
younger than 1 Myr. 
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Figure 22. Theoretical isochrones in the range 1-100 Myr for the labelled chemical compositions with a - 1 .68 and Xd =4 10 
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Table 3. Summary of the models and isochrones available in the database. The table lists the initial deuterium abundance Xd, the 
mixing length parameter (a), the primordial abundance of helium (Yp), the helium to metal enrichment ratio (AY/AZ), the initial 
helium (Y), and metal (Z) abundance. 
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Solar chemical composition Y = 0.2533 Z = 0.01377 a = 1.68 Xd = 2.0 • 10-= 



chemical composition, we have made available 43 evolutionary 
tracks computed with three different values of the mixing-length 
parameter a - 1.2, 1.68, and 1.9, in the mass range 0.2 - 7 M©, 
and 36 isochrones, in the range 1 - 100 Myr. For Z > 0.008, two 
different initial abundances of deuterium have been adopted. The 
database is available on the wetQ Models are compared in de- 
tail with other computations available in the literature, while the 
comparison with selected observational data shows good agree- 
ment. 
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